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We present a statistic for the detection of stochastic gravitational wave backgrounds (SGWBs) 
using radiometry with a network of multiple baselines. We also quantitatively compare the sensitiv- 
ities of existing baselines and their network to SGWBs. We assess how the measurement accuracy 
of signal parameters, e.g., the sky position of a localized source, can improve when using a network 
of baselines, as compared to any of the single participating baselines. The search statistic itself is 
derived from the likelihood ratio of the cross correlation of the data across all possible baselines in 
a detector network and is optimal in Gaussian noise. Specifically, it is the likelihood ratio max- 
imized over the strength of the SGWB and is called the maximized-likelihood ratio (MLR). One 
of the main advantages of using the MLR over past search strategies for inferring the presence or 
absence of a signal is that the former does not require the deconvolution of the cross correlation 
statistic. Therefore, it does not suffer from errors inherent to the deconvolution procedure and 
is especially useful for detecting weak sources. In the limit of a single baseline, it reduces to the 
detection statistic studied by Ballmer [ ] and Mitra et al. [ ]. Unlike past studies, here the MLR 
statistic enables us to compare quantitatively the performances of a variety of baselines searching 
for a SGWB signal in (simulated) data. Although we use simulated noise and SGWB signals for 
making these comparisons, our method can be straightforwardly applied on real data. 

PACS numbers: 95.55.Ym, 04.30.Db, 98.80.-k, 97.80.-d 



I. INTRODUCTION 

Just like the discoveries of the cosmic microwave back- 
ground and pulsars in the electromagnetic spectrum, a 
discovery of unknown sources by Earth-based detectors 
such as LIGO and Virgo in the gravitational wave (GW) 
spectrum by serendipity is an interesting prospect. The 
LIGO Scientific Collaboration and the Virgo Collabora- 
tion are addressing it by searching for both transient sig- 
nals, or "bursts," and long-duration signals in the data 
from their detectors. Here, we focus on a subset of the 
latter type that can be modeled as a stochastic back- 
ground. The search for an isotropic stochastic GW back- 
ground has caught significant attention due to its cos- 
mological significance. This primordial GW background 
is a direct probe of cosmological inflation [ ]. However, 
the astrophysical background, arising in the nearby Uni- 
verse [5], e.g., from an unresolved superposition of GW 
signals from multiple sources, such as low-mass x-ray bi- 
naries or, even, coalescing compact objects, is possibly 
much stronger than the primordial background and is 
anisotropic. 

A variety of data analysis strategies to search for an 
anisotropic GW background have been proposed and im- 
plemented in the past [6-11]. These searches are usu- 
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ally performed in two types of bases in the sky, namely, 
the pixel and spherical harmonic bases. Use of the ra- 
diometer technique for searching a GW background was 
proposed in Ref. [ 2] and was implemented in the pixel 
basis on data from LIGO's fourth science run [ ]. An 
elaborate study of this method, including the maximum- 
likelihood (ML) estimation of the true anisotropy of GW 
background by deconvolving the observed sky map, was 
presented in Mitra et al. [ ]. Even though the pixel-based 
search is promising and simpler to understand, it is not 
the best basis for probing sources with angular spreads 
greater than the angular resolution of the GW radiome- 
ter. The spherical harmonic basis is expected to be bet- 
ter suited for detecting such sources [ ]. Past attempts 
at probing the GW anisotropy in the spherical harmonic 
basis were essentially studies of the periodic modulation 
of the observed background in the detector baselines. Re- 
cently, a general ML formalism was developed to search 
for the GW anisotropics in any basis, including the spher- 
ical harmonic basis, using a network of ground-based GW 
interferometers [ ]. The pixel-based search is a specific 
application of this formalism. One of the main goals of 
this paper is to perform a thorough comparison of the 
expected performances of individual baselines and the 
whole network in detecting an astrophysical stochastic 
gravitational wave background (SGWB) and in estimat- 
ing its parameters. The pixel basis is used for this study. 

Even though a pixel-based search is optimal for a lo- 
calized source, the resolution of the source is limited by 
the length of the radiometer baselines, the orientation of 
the detectors, and their individual sensitivities. Prob- 
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ing a stochastic GW background with energy distributed 
across the pixelated sky demands a statistically mean- 
ingful integration of the energies received in every pixel. 
In order to accomplish this, we extend the maximized- 
likelihood ratio (MLR) statistic for a single baseline to 
incorporate a network of detectors or, equivalent ly, mul- 
tiple baselines. The rest of the paper is devoted to study- 
ing the performance of individual GW detector baselines 
and the whole network by comparing different figures of 
merit for their performance, e.g., sensitivity, accuracy in 
localizing sources, sky coverage, and faithful extraction 
from the data of the sky distribution of a stochastic back- 
ground. 

The paper is organized as follows: In Sec. II, we de- 
velop and study the efficiency of an optimal all-sky search 
statistic for anisotropic SGWBs that obviates the solving 
of the inverse problem, which may not always be well- 
posed. In Sec. Ill, we compare the performance of a net- 
work with that of its individual baselines using a variety 
of figures of merit. In Sec. IV, we conclude by summa- 
rizing the implications of this work on ongoing SGWB 
searches and by highlighting future directions in GW ra- 
diometry. 



II. OPTIMAL SEARCH STATISTIC 

A. Statistical properties of the signal and detector 

noise 

In the transverse traceless gauge, the spatial part of 
the metric perturbations due to a SGWB can be written 
as a superposition of plane waves 
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dCl e*(£l) h A (f,Cl)e i2 ^ t+fl ^ , 

(2.1) 

where a and b are spatial indices, e^ b (Q) are the com- 
ponents of the gravitational wave polarization tensors, 
and Cl is a unit vector pointing in the direction of 
wave propagation. Note that the Einstein summation 
convention for repeated indices has been applied to the 
polarization index A = {+, x}. Since the GW strain 
h a b(t, r) is real, the complex Fourier amplitudes fiA{f, ^) 
satisfy the reality condition, 1ia(— f, = h* A (f,Q). 

The polarization tensors can be defined in terms of the 
spherical polar coordinates, G [0,7r] and <p G [0, 2tt), on 
the sky. Let as follows 

Cl = sin cos (j) x + sin sin <\> y + cos z , 
m = cos cos (j) x + cos 6 sin <j> y — sin z , 
n = — sm(f)x + cos (j)y , (2.2) 

such that {m, n, ft} form a right-handed system of unit 
vectors. The axes are defined as follows: For a fixed but 
arbitrarily chosen origin of time t = 0, x is directed to- 
ward the intersection of the equator and the longitude 



(j) = 0, z points at the Celestial North Pole, and y is 
chosen orthogonal to the previous two axes, forming a 
right-handed triad. Then, the polarization tensor com- 
ponents e£ b (U) are defined as 



e ab&) = m a rhb ~ n a n b , 
e ab&) = mafib + h a m h , 



(2.3) 



in the aforementioned right-handed orthogonal basis. 

Understanding the signal excited in an interferometric 
detector by a SGWB is helped by specifying the detec- 
tor's location and orientation in the above orthogonal 
basis. Let the I^ 1 GW detector be located at 17 (t), and 
let Xi(t) and Yi(t) be the unit vectors pointing along 
its arms. These three detector location and orientation 
vectors are all time-dependent due to the Earth's rota- 
tion. Then, the components of the 1^ detector tensor 
are given by 



df{t) = - xf(t) x b j(t) - Y?(t) y/(t) 



and 



h I (t) = h ab (t,r I (t))df(t), 



(2.4) 



(2.5) 



is the strain in it due to the SGWB. 

The response of a detector to the polarization compo- 
nent A of a wave incident from direction Cl is given by 
the antenna-pattern function 



(2.6) 



where we assumed the Einstein summation convention 
over the repeated indices a and b. Contracting (2.1) with 

the detector tensor the GW strain signal in the I 
detector can be expressed as 
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df / dnh A (f,n)F?(n,t)e 
Js 2 



Virf(t+£lTi(t)/c) 

(2.7) 

in terms of the antenna-pattern functions. 

The Fourier components of the strain ^) describ- 

ing a stochastic GW background are random variables 
whose expectation values define the statistical properties 
of the background. Without loss of generality we assume 
that these components have zero mean: 



(h A (f,n)) = o, 



(2.8) 



where the angular brackets denote statistical average. In 
the presence of a signal, the time series of the 1^ detec- 
tor's output xi(t) is a sum of the GW signal hi(t) and 
the detector noise nj(£): 



xi{t) = h T (t) +7l/(t). 



(2.9) 



Statistically, the gravitational wave strain hj(t) are un- 
corrected with the detector noise, implying that 



(h I {t)n J {t')) = 0. 



(2.10) 
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We also assume that the noise is Gaussian with zero 
mean, i.e., (n/(£)) = 0, and is uncorrelated in different 
detectors, namely, 



(ni(f)nj(f')) = -d(f-f')S IJ ^ I) (f): 



(2.11) 



where is the one-sided noise power spectral density 
(PSD) of the 7 th detector. The last assumption is not 
unreasonable when the detectors are widely separated 
across the globe. 



where the beam function B A k , is analogous to the point- 
spread function that maps the power in the object (or 
sky) plane to that in the image plane. Above n k is the 
noise in the k th sky position, and S k is termed as the 
dirty map [ ]. We define as a vector, with V k A ^ as 

its k th component, and IS a as a matrix, with B Ak , as its 
(k, fc') th element. 

C. Detection statistic 



B. Cross correlation statistic 

Since the targeted source is stochastic, we search for 
its GW signal by looking for correlated patterns in the 
data of two or more detectors after accounting for time 
delays and detector responses consistent with a given sky 
location. This is done by cross-correlating the data xi(t) 
from the detectors, taken in pairs, with a sky-position- 
dependent time- frequency filter Q k (t;f), labeled by the 
sky-position index k. The cross correlation statistic com- 
bined for the observation period T for the data xi^(t) 
from two detectors or, equivalent ly, for a baseline is de- 
fined as 



S k = 4 At 



t=0 J ~< 



df xUt;f)x 2 (t;f)Q k (t;f), (2.12) 



where xj(t]f) is the short-term Fourier transform of 
xi(t), over time interval At, and is defined as in Ref. 



pt+At/2 

xj{t-f) := / dt'xjit'je- 2 ^' . (2.13) 

Jt-At/2 

The filter that maximizes the signal-to-noise ratio 
(SNR) associated with this statistic is a scalar, square- 
integrable function on the sky [ ] and, hence, can be 
resolved linearly in an appropriate basis, such as a pixel 
basis or the spherical harmonic basis. In the former 
case, k is the pixel index. 

Let an astrophysical GW background be modeled such 
that the Fourier components of its GW strain h A (f) of 
polarization A from the k th sky-position obey 

{h k A *(fM,(f')) = 5aa> S(f - f')5 kk 'V k {A) H(f) , (2.14) 

where V k A \ is a dimensionless measure of the signal 
strength, and H(f) is its two-sided power spectral den- 
sity, with units of Hz -1 [2]. Here, we assume the signal 
to be stochastic and uncorrelated in the two polariza- 
tions, different frequencies, and different sky locations. 
In the presence of a signal in the detector data, the cross 
correlation statistic is 



S k =B k +k ,V^ )+ B k k ,V^ )+ n k 



To get a single detection statistic, one must combine 
the measurements of S k for all k. When the detector 
noises are Gaussian and uncorrelated, an assumption 
borne out to sufficient approximation for our purposes, 
the n k are Gaussian with a nontrivial covariance matrix, 
N, determined by the beam functions. The exact form 
of N is discussed below. 

If an astrophysical GW background signal, character- 
ized by the pixel-strength vector P, is present in the data, 
then the probability density function of the radiometer 
output S is given by 

p(S\V) = (2ir)- N ^ /2 

x exp[-i((S-B-P) T -N- 1 -(S-B-P) 

+ Tr[lnN])] , (2.16) 

whereas in the absence of a signal it is 

p(n) = (27r)- A/ "p ix/2 exp[-i(n T - INT 1 • n + Tr[lnN])] . 

2 (2.17) 
By the Neyman-Pearson criterion, the optimal detection 
statistic is the likelihood ratio p(S\V)/p(n) [14]. 

For an unpolarized background from a source dis- 
tributed across multiple pixels and quantified by the 
signal-strength vector V = P(+) = ^(x)> the log- 
likelihood ratio maximized over V = IIP II is 



k' 



X 



S k (N- 1 ) kk ,(B.V) 

s k v k 



V^B qr V r 



(2.18) 



where TP is the unit vector along The beam matrix 
for an unpolarized source is given by 

(2.19) 



B pq — B+p q + B 



Xpq i 
T 



= SAf At^2 i)T(Cl p , t) 



t=o 



fu 



e 27rif(Q q -Q p )-Ax(t)/cQ^j^ 

f=fi 



(2.15) 



(2.20) 
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where T(Cl q ^t) is the time-varying baseline antenna pat- 
tern, and G(t, f) is a measure of the spectral strength of 
the source relative to the baseline's noise PSDs: 



T(Cl,t) := F+(£l,t)F+(fl,t) +F 1 x (^,t)F 2 x (^,t), 



G(tJ) :-- 



f(l)(*,/)f(2) (*,/)" 



(2.21) 



In the weak-signal limit the noise-covariance matrix is 
approximately equal to the beam matrix, 



(2.22) 



Its diagonal elements inform us about the sensitivity of 
the network to the different pixels in the sky for an 
SGWB with PSD H(f). 

The statistic A is the maximized (log-) likelihood ra- 
tio for a single-baseline SGWB search and is the same 
statistic introduced in Appendix C of Ref. [ ]. Here, it 
is expressed specifically in terms of quantities defined in 
the pixel basis. It has zero mean and unit variance in 
the absence of a signal. When a signal is present in the 
data and its parameters are matched exactly by the tem- 
plate's, the mean of the statistic is 



(2.23) 



The variance of the statistic remains unchanged. One 
can extend this single-baseline statistic to the case of a 
multibaseline network. That statistic arises directly from 
maximizing the log-likelihood ratio for a network and is 
given by 



k' 



ExUs x k v k 



(2.24) 



where X is the baseline index and the subscript M high- 
lights that this MLR statistic is for a network of baselines. 

The MLR statistic is a detection statistic for SGWBs in 
the same manner as the standard matched-filter statis- 
tic is for deterministic GW sources. The latter is also 
obtained by maximizing the likelihood ratio with respect 
to the strength of the deterministic source. Searching for 
a signal from a deterministic source involves maximiz- 
ing the matched-filter statistic over a bank of templates 
defined on the signal's parameter space. For SGWBs, 
as well, the detection statistic can be the MLR, max- 
imized further with respect to different SGWB models 
given by *P, perhaps parametrized by a smaller number 
of parameters than the number of components of TP. This 
is in contrast to the existing searches for anisotropic GW 
backgrounds. Past dirty- map-based searches precluded 
the presence of a signal by demonstrating that the map 
is consistent with a Gaussian distribution, up to statis- 
tical fluctuations allowed by the number of independent 



"samples" on the sky [ ]. However, they did not provide 
a confidence level for the presence or absence of a broad- 
band or spatially extended signal. A better approach 
is to solve the inverse problem in an orthogonal basis, 
namely, the pixel [2] or spherical harmonic [3] basis. This 
yields an estimate of the background, i.e., a "clean" (de- 
convolved) map and the corresponding noise-covariance 
matrix. However, this approach depends heavily on how 
well-posed the inverse problem is and how accurately it 
can be implemented numerically. Consequently, a detec- 
tion statistic constructed on the deconvolved data can be 
affected by similar maladies. To work well, the inverse 
problem requires that the network of interferometers is 
sufficiently nondegenerate, which is not always the case. 
Indeed, the deconvolution procedure can enhance spatial 
noise correlations and, sometimes, even introduce arti- 
facts, thereby adversely affecting parameter estimation 
and signal detection by such a procedure. 

As we prove here, the detection problem does not re- 
quire a well-posed inverse problem and exists even for 
a degenerate network. A detection statistic is best de- 
fined on the dirty map, as opposed to the clean map. 
As an added advantage, a dirty-map-based statistic is 
faster to compute, since it obviates the computational 
overhead required for obtaining the clean map. While 
it is possible to use an arbitrary sky model, such as the 
one- dimensional basis P, and estimate the strength of 
the SGWB for that particular model, the MLR statistic 
in Eq. (2.18) provides a well-understood construct that 
can be maximized over a set of parameters for selecting 
the model that best fits the data. 

To elaborate further on the way the new statistic 
works, let us consider the example of a directed pixel- 
space search, which is performed for only one source and 
assumes that the angular extent of the source is, at most, 
one pixel. In the standard radiometer search [ , ], the 
dirty map S k is computed for each pixel k in the sky. 
The (signal part of the) dirty map is generally peaked 
at the source pixel and has broad structures, including 
large negative patches, around it. One way of inferring 
the presence or absence of a source in this image requires 
deconvolving it. However, as we show later, deconvolu- 
tion of a relatively weak source can result in a clean map 
with significant errors, especially when the sky is divided 
into around 3000 pixels or more. (A network resolution 
of several square degrees requires a few thousand pix- 
els across the sky.) Also, computing the noise-covariance 
matrix can be numerically challenging. The MLR is a 
good choice in this situation, since it combines all the 
pixel values to provide a single number for the detection 
statistic that is simple to use in drawing inferences on 
the presence or absence of a signal in the network data. 

If a parametrized model of the background is avail- 
able, one can construct the likelihood-ratio statistic from 
the dirty map and maximize it over the parameter space. 
The maximized likelihood-ratio statistic can also be used 
to perform a more advanced blind search, where no prior 
information is available about angular distribution of the 
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power in the SGWB. For each basis component, one can 
assert that only that basis component is present in the 
signal and compute the statistic with the correspond- 
ing sky model. Thus, using the dirty maps of "point 
estimate" or SNR obtained by the existing radiometer 
search, our prescription takes one step forward and can 
provide a map of likelihood ratios, which is statistically 
a more robust and meaningful quantity, given a set of 
highly correlated observations. 

The construction of the MLR statistic on a dirty map is 
simple. Equation (2.18) shows that it is the scalar prod- 
uct of the observed map S k and a sky-model-dependent 
normalized "template." The template is proportional to 
B • P, which is the expected signal in the dirty map for 
a sky model TP. The inverse of the noise-covariance ma- 
trix is the metric in the pixel space. The sky model can 
be defined in a straightforward way. For instance, to 
search for a point source localized to a single pixel, one 
would use a TP with all but one component, namely, the 
component corresponding to that pixel, set to zero. In- 
deed, B • V is simply the point-spread function of the 
pixel with the nonzero component of TP. Also notice that 
the inverse of the noise-covariance matrix, being propor- 
tional to the beam matrix, cancels out algebraically in 
the expression for the MLR. Therefore, unlike for decon- 
volution, for MLR construction the computation of this 
matrix is not needed. Otherwise, the latter procedure 
would have been computationally similar to solving the 
inverse problem, avoiding which is one of the main moti- 
vations for this work. 

To complete the discussion, we note that the construc- 
tion of a MLR statistic is not limited to dirty maps and 
can be implemented for clean maps. A clean map can be 
expressed as 

V = V + n c , (2.25) 

where V is the true sky map, and n c is Gaussian 
noise with covariance XI, which is related to the dirty- 
map noise-covariance matrix through the relation X) = 
(B T N- 1 B)- 1 . Therefore, following the same procedure 
as that for the dirty map, one can write the MLR statistic 
for a clean map as 

V • E" 1 • V 

A c = , (2.26) 

W • X) -1 • V 

and, thereby, obtain model-based or blind likelihood- 
ratio maps. 

D. Performance of optimal detection statistic 

We numerically study the performance of the optimal 
statistic and compare with the existing method. We use 
the LIGO 4km detectors located in Hanford (HI) and 
Livingston (LI). Unless otherwise stated, the noise PSDs 
of all detectors are taken to be their (smoothed) first- 
generation design sensitivities [28, 29]. The frequency 



band considered here spans 40 — 1024Hz, with a bin size 
A/ = 1Hz. The source PSD is taken to be a constant, 
H(f) = 1.516 x 10" 48 /Hz. Note that the spectral index 
of the source PSD has a significant effect on the resolu- 
tion of the network. Predictions from astrophysical and 
cosmological models suggest the nominal range of the 
spectral index to be between —3 and 1. The higher the 
spectral index, the higher the resolution, and the more 
computationally expensive the directed search. We take 
the spectral index to be zero here by setting H(f) as a 
constant. 

The directed search is performed by dividing the (sim- 
ulated) strain data from all detectors into segments with 
a duration of 192 sec. The noise is taken to be station- 
ary. The sky is tessellated into 3072 pixels by using the 
Hierarchical Equal Area isoLatitude Pixelization [15, 16]. 
The choice of the signal integration duration is taken to 
be a sidereal day, which leads to the azimuthal symme- 
try of the baseline sensitivities and sky resolutions. The 
justification for choosing the above parameter values can 
be found, e.g., in Mitra et al. [ ]. 

We first construct simulated data sets of two kinds, one 
with only noise and the other with a weak signal from a 
"polar-cap" source added to that noise. The sky map of 
the latter case is shown in the first plot in Fig. 1. We 
make dirty maps for these two cases using single-pixel 
source templates for each of the 3072 pixels. These maps 
are essentially the maps of SNR for the directed search, 
as can be seen by substituting 

fik = f 1 for target pixel , ^ ^ 

1 for remaining pixels , \ • J 

into Eq. (2.18) [2]. The dirty maps for both cases look 
very similar, and only one of them, namely, the one for 
the polar-cap source, is shown in the second plot of Fig. 1. 
Not surprisingly, it is also similar to the dirty map pre- 
sented in Abbott et al. [ ] for real data. Indeed, the 
MLR values over the dirty-map pixels appear to follow a 
normal distribution, as shown in Fig. 2. This is consis- 
tent with the distribution presented in Abbott et al. [13]. 
Following that reference, we also plot the la error enve- 
lope around the Gaussian fit for 400 degrees of freedom 
and observe that the tops of every bar in the histogram 
lie within that envelope. The important point to note 
here is that, for both noise-only and weak-signal (polar- 
cap) data sets the distributions in Fig. 2 are very similar 
and consistent with a normal distribution. 

If we now pretend that we know the broad shape of the 
GWB sky and use the V k of the polar-cap signal as our 
template for computing the MLR statistic, we find that 
the above two cases can be distinguished better: In the 
noise-only case, the MLR statistic equals —0.023, while 
with the weak polar-cap signal it is 1.400, which is sig- 
nificantly larger than the former, as explained below. To 
corroborate this claim, we computed the MLRs on an 
ensemble of 4000 realizations of noise, with and without 
the weak polar-cap signal. In Fig. 3, we show the dis- 
tribution of the MLR statistic (A) for noise-only (top) 
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(c) Clean map 

FIG. 1: The V k sky map of a weak polar-cap GWB source 
is shown in (a). The dirty (b) and clean (c) maps for this 
source were constructed using the radiometer algorithm for 
the LIGO H1L1 baseline. The last two maps for this weak 
source are visually very similar to those for the noise-only 
case (which is not shown here). 



and weak-signal (bottom) cases. Clearly, the noise-only 
A values are normally distributed with a zero mean, and 
the weak-signal A values are normally distributed with a 
mean of ~ 1.4<j, where a ~ 1. This experiment confirms 
that, given our assumptions on the detector noise and the 
signal, the MLR statistic can considerably enhance the 
detectability of a weak diffuse stochastic background, if 
a reasonable model of the background is available. This 
is true even when the distribution of the dirty-map pixel 
values for that source is close to zero-mean Gaussian. 



III. PERFORMANCE OF MULTIBASELINE 
RADIOMETERS 
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(a) Noise 
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(b) CAP 

FIG. 2: Histograms of the dirty maps for the noise-only 
data set (top) and weak polar-cap signal (bottom), which is 
depicted in Fig. 1, are shown here. These two histograms 
are consistent with that of zero-mean Gaussian data (solid 
curved line), up to la errors: Following Abbott et al. [13], 
the la error boundaries for 400 degrees of freedom have also 
been overlaid for a consistency check. 



In this section we define a set of figures of merit to 
compare the performance of a network of baselines with 
that of its individual baselines for a directed search of a 
SGWB. 

A single figure of merit may not suffice in capturing 
all the attributes of a baseline or a network of detec- 
tors benefiting a SGWB search. A certain baseline or 
network configuration can have good sensitivity if the 
detectors are optimally orientated, but have poor resolu- 
tion if they are proximally located. Optimally oriented 
detectors may be very sensitive to certain anisotropy 
modes, but insensitive to others, making the estimation 
problem highly degenerate. On the other hand, a net- 
work of detectors that are oriented differently may have 



moderate, yet uniform, sensitivity to all spherical har- 
monic modes of a SGWB, thereby mitigating the ill- 
posedness of the estimation problem. Such a network, 
however, will perform worse than one where all the de- 
tectors are aligned similarly in a low- frequency, all-sky 
isotropic search. Therefore, the relevance of a figure of 
merit is determined by the kind of search one is under- 
taking. Here, we propose a set of figures of merit that 
are relevant to current searches of anisotropic stochastic 
background and that are special cases of the general ML 
framework presented in Thrane et al. [ ]. 

For the numerical simulation studies below, we use 
the same detector characteristics as mentioned in Sec- 
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(a) Noise 
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(b) CAP 

FIG. 3: Monte Carlo simulations with 4000 noise realiza- 
tions were performed to study the performance of the MLR 
statistic. The distribution of the MLR statistic obtained for 
dirty maps generated from noise-only (top) and weak polar- 
cap injection (bottom) are shown. Clearly, the MLR statistic 
detects the signal at ~ lAa level. 



tion II D, but we now include the Virgo detector (VI) 
in Cascina, Italy to construct a three-baseline network. 
The baselines and their network are named by concate- 
nating the symbols for the participating detectors; e.g., 
the Hanford-Livingston baseline is termed as H1L1, and 
the network of the above three detectors is termed as 
H1L1V1. 



A. Sensitivity 

The first figure of merit is the "sensitivity" of a net- 
work and is motivated by a similar quantity defined in 
Cella et al. [ ] for the all-sky isotropic search. In prac- 



tice, a greater sensitivity implies a better confidence level, 
at which detection can be made or upper limits can be 
inferred. 

We define the single-baseline sensitivity for a directed 
search as the expectation value of the MLR in Eq. (2.23) 
for a SGWB source with V set to unity, 

Sensitivity = \j (B • V) k (N~ 1 )kk , (B • V) k ' , 

= \lv k B kk >V k ' . (3.1) 

The sensitivity can be expressed in the spherical har- 
monic basis as follows: 

Sensitivity = ^V lrn B lmVm 'V V m> , (3.2) 

where 

V lm = [ dClV{h)YC m (Sl) , (3.3) 



Blml'rn' = J J dQdCl' Y lm {Cl)B{(l,Cl')Y Vm ,{(l') . 

(3.4) 

Owing to the statistical independence of the baselines, 
the multibaseline sensitivity squared is the sum of squares 
of the individual baseline sensitivities, as was also noted 
for the isotropic-background baseline sensitivities in Cella 
et al. [17]: 

Sensitivity^ = ^ Sensitivity^ . (3.5) 
x 

For an unpolarized background from a single pixel, say, 
labeled fc, and with V r = 8 r ^ k \ the sensitivity expression 
simplifies to 

Sensitivity (/c) = ^(^(N" 1 )^^*) = y/BWW . 

(3.6) 

Unless otherwise mentioned, there is no sum over the 
repeated parenthetic indices in this paper. 

In the top panel of Fig. 4, we compare the sensitivi- 
ties of the baselines and the whole network as a function 
of declination. (As noted above, the sensitivities are az- 
imuthally symmetric.) For a fair comparison, we also 
replot them after weighting them with the cosine of the 
latitude, in effect, to assign equal weight to every pixel 
on the sky. It is clear that the H1L1 baseline has much 
better sensitivity due to the similar orientations of the 
two detectors. Still, inclusion of Virgo, which is oriented 
quite differently relative to HI and LI, improves the sen- 
sitivities of the network by ~ 10% (which corresponds 
to an increase in the observational volume by ~ 30%), 
especially in the regions where the H1L1 baseline does 
not perform well. However, this network improvement 
is highly superseded by all other performance improve- 
ments indicated by corresponding figures of merit intro- 
duced in this section. 
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In Fig. 5 we plot the narrowband (5Hz) sensitivities 
at two locations, namely, the Celestial North Pole (top) 
and the equator (bottom). Performance improvement of 
a network for a narrow band search at high frequencies is 
better than the (frequency integrated) broadband search. 
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FIG. 4: The sensitivities (top) and their area- weighted coun- 
terparts (bottom) of three different baselines and their net- 
work are plotted as functions of the declination of a single- 
pixel SGWB source. The source PSD (H(f) = 1.516 x 
10 _48 strain/Hz) is chosen such that it has maximum SNR= 
10 in the H1L1 baseline. (Note that the source parameter V 
is set to unity for these plots.) The signal band considered 
here is 40-1024 Hz. 



B. Sky coverage 

In a directional search, the main advantage of a net- 
work lies in the fact that it vastly improves sky coverage, 
which, in turn, leads to better parameter estimation, in- 
cluding localization and sky-map reconstruction consid- 
ered later in this section. In this subsection, we illustrate 
the advantage of using a network of detectors, as com- 
pared to using its individual baselines, to this end. 

In general, one radiometer baseline cannot sample the 
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FIG. 5: The sensitivities of three single baselines and their 
multibaseline network plotted as functions of the central fre- 
quency of the source band. The source is chosen to have a 
constant H(f) — 1.516 x 10 -48 strain/Hz and a band width of 
5Hz. The top panel represents the sensitivities at the celestial 
poles, and the bottom panel represents those at the celestial 
equator. 



whole sky uniformly; the measurement errors in some 
parts of the sky are much worse than those in the other 
parts. Introduction of new baselines with different ori- 
entations improves filling in these "holes" by scanning 
the sky with different antenna-pattern functions. In the 
first three plots in Fig. 6, we show the standard devia- 
tion in measuring the dirty map by the three individual 
baselines, namely, H1L1, HI VI, LI VI, and their network 
H1L1V1 [ ]. The azimuthal symmetry mentioned before 
is explicitly observed here. The H1L1 baseline has the 
least deviation at most declinations, due to the optimal 
orientations of the HI and LI detectors. Again, H1V1 
and LI VI baselines have low deviation in the regions 
where H1L1 does not perform well. Since the dirty maps 
from different baselines are operationally combined with 
an inverse- noise- variance weight, the harmonic mean of 
the variances provides the effective variance of the com- 
bined dirty map. The last plot in Fig. 6 shows the ef- 
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fective deviation for a network of detectors. Clearly, the 
deviation now has smaller spread and also, by construc- 
tion (harmonic mean), the deviations are smaller than 
those for the individual baselines. 




0.105 0.11 0.115 0.12 0.125 0.13 0.25 0.3 0.35 



(a) H1L1 (b) L1V1 



numerical errors in the anisotropy estimation problem 
at "high" resolution, i.e., near or beyond the diffraction 
limited resolution, and significantly improves source lo- 
calization accuracy. 

Singular value decomposition of the Fisher informa- 
tion matrices provides a more quantitative verification of 
the above claim. Figure 8 shows the singular values of 
the Fisher matrices for the individual baselines and the 
whole network. The LIGO baseline has very small sin- 
gular values at higher resolutions (dashed curved line), 
which implies that estimation of anisotropy at those res- 
olutions is an ill-defined problem. The network reduces 
the difference between high and low singular values and 
regularizes the inverse problem at high resolution (solid 
curved line). 
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(c) H1V1 (d) H1L1V1 

FIG. 6: The standard deviation of dirty maps measured by 
the three LIGO- Virgo baselines and the full network for a 
constant H(f) are plotted here. The aim of this figure is 
to show how individual baselines complement each other, 
thereby making the "effective deviation" (i.e., the square root 
of the harmonic mean of variances of the individual baselines) 
of the combined map observed by the network nearly uniform. 
Note that the color scale in the network plot has lesser spread 
than the individual baselines. The absolute scale of the maps 
depends on the normalization of the filter, and only the rel- 
ative scale is important here. The azimuthal symmetry is 
present because we are considering a whole sidereal day's ob- 
servation, with stationary noise. 



Most importantly, a network also complements the 
single-baseline observations in terms of angular resolu- 
tion. The beam functions for each radiometer baseline 
are highly asymmetric, which means that a given posi- 
tion on the sky is probed with quite different angular 
resolutions in the tangential directions. To illustrate this 
aspect, the typical beam functions for the three LIGO- 
Virgo baselines in the direction of the Virgo cluster are 
shown in Fig. 7. If we consider the beam for the H1L1 
baseline, the sensitive part of the beam is similar to a 
highly eccentric ellipse, suggesting that the angular reso- 
lution along the minor axis is much finer than that along 
the major axis. The beams for the baselines in a network 
involving the Virgo detector provide better resolution due 
to the longer baselines: The beams are finer along the 
major axis of the H1L1 beam, thus complementing the 
H1L1 observation, which is a major motivation for using 
a network. This, in turn, improves the condition number 
of the Fisher information matrix, thereby, reducing the 



C. Parameter accuracy 

An important figure of merit for a directed search is 
how well a point source can be localized or its other pa- 
rameters be constrained. In a noise-dominated mapping 
experiment, it is not easy to identify sources in the ob- 
served images. If there are candidate sources that have 
been modeled by other astronomical observations, one 
can utilize that information to detect or constrain pa- 
rameters of such sources. If the source was very accu- 
rately modeled, the optimal strategy would be to design 
a specific search focused on that source. But, in practice, 
with very limited knowledge of sources, the optimal strat- 
egy would be to vary the parameters within a reasonable 
range and maximize the log-likelihood ratio. 

One of the main advantages of the MLR statistic is that 
it allows estimation of parameters of the source, given 
a model. For example, if there is an extended source, 
such as a cluster of galaxies, with an angular scale com- 
parable to the resolution of the radiometer, and there 
exists a reasonable model for its mass distribution, one 
can maximize the log-likelihood ratio to find the center 
of the cluster [33]. Even for a blind search, this method 
may prove to be advantageous to perform a finer search 
around the poorly estimated parameters of a potential 
candidate source. 

In this section, we assess the accuracy with which a 
pointlike (single-pixel) source can be located using a net- 
work of GW detectors as compared to its individual base- 
lines. The parameter estimation accuracy is deduced 
from the elements of the Fisher information matrix [14]. 

For an unpolarized background from a single pixel, la- 
beled /c, and with V r — S r ^ k \ the single-baseline detec- 
tion statistic follows from Eq. (2.18) to be 

S^N- 1 )^^ 
\k) = / > (3-7) 

ytef fe) (N-i) ra ^ fe) 

which can be interpreted as the inner product of the data, 
S, and a unit-norm template Hence, the match [18] 
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FIG. 7: The beam functions for the three LIGO- Virgo baselines for H(f) = constant are shown here. This figure illustrates 
that different baselines also complement each other in terms of angular resolution along different tangential sky directions. 



Singular Values of the Fisher Matrices of SGWB Directed Search 
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FIG. 8: Singular values of the Fisher matrices for individual 
baselines and the whole network are plotted in this figure. 
The LIGO H1L1 baseline (dashed curved line) has very small 
singular values at high resolutions; therefore, estimation of 
anisotropy at those resolutions is an ill-defined problem. The 
network (solid curved line) makes the singular values much 
more uniform, thereby regularizing the inverse problem. 



between the unit-norm templates for the k th and the k' th 
pixels is 



M 



(*)(*') 



(3.8) 



where the inner products are all defined in terms of N _1 . 
Define &(k) = {fik^k} as the pixel coordinates, where 
ji k = cos#/e, with Ok and fa being the declination and 
right ascension of the k th pixel, respectively. Since the 
match has a maximum value of unity at k' = fc, one 
can expand M in a Taylor series about A/i^) = and 



M « 



as 
1 



1 



d 2 M 



2 yaefaaefa 



© 



(k'y 



=0 



(fc) 



(SNR)* AO W AO W = 



(3.9) 



where T 



(k)fiu 



r /ii /(0(/ c )) are the components of the 
Fisher information matrix T^), and (SNR)^) is the 
signal-to-noise ratio in the k th pixel. 

For large SNR, the error variance-covariance matrix 
obeys 



Aef fc) Ae^ fe) 



(3.10) 



The estimation error in the measurement of the sky- 
position solid angle (in steradians) is given by [ ] 

27r^/<(Acos0 (fe) )2) ((A0 (fc) )2) - (Acos0 (fe) A0 (fe) ) 2 . 

(3.11) 

The Fisher information matrix for multiple baselines is 
just the sum of the Fisher matrices for the individual 
baselines, 



[r 



X(k)jj,u ; 



where X is the baseline index, and T x 



(k)fiu 



(3.12) 



is the 



Fisher information matrix of the X th baseline, as given 
in Eq. (3.9). Hence, the error variance-covariance matrix 
for the network is 



1\ V 



AO^AOf, 



(fe)^(fe)/ 



(3.13) 



for large SNR. Therefore, the 1<j estimation error in solid 
angle for locating a pixel source with the multibaseline 
network is expressed as 
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[An 



2^^((Acos^ ) ) 2 ) Ar ((A(/> (fe) )2) Ar -(Acos^ ) A(/) (fe) ) 



2 

AT 



(3.14) 



r 



Note that this error diminishes with SNR as 1/SNR 2 , 
i.e., localization is more accurate at higher SNR. 

We present the source-localization errors for the indi- 
vidual LIGO- Virgo baselines and the network in the top 
panel of Fig. 9. We also show the corresponding area- 
weighted plots obtained by multiplying these errors with 
the cosine of the latitudinal angle in the bottom panel of 
Fig. 9. The network clearly outperforms individual base- 
lines by about 1 order of magnitude or more for almost 
all declination angles. 
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FIG. 9: The la error (top) and area- weighted la error (bot- 
tom) in the solid angle for locating a source with three single 
baselines and the whole network. The network accuracy is 
better by about 1 order of magnitude or more at most of 
the declinations. Note that an error of lsr ~ 3282.80635 sq- 
degrees, and that the error here decreases as 1/SNR 2 . 



The primary focus of this analysis was to obtain statis- 
tics that are based on dirty-map constructs. However, it 
is straightforward to extend it to be applicable to clean 



maps. Similarly, although we considered broadband sig- 
nals, it is possible to easily extend our study to narrow 
band signals. 



D. Map making 

Finally, we compare the quality of sky maps made 
by the individual baselines and their network since they 
are among the primary products of anisotropic searches. 
Here, we consider two figures of merit, namely, the MLR 
statistic and the normalized mean square error (NMSE) 
for comparing maps. To compare the dirty maps, we 
use the MLR statistic introduced earlier, and, to com- 
pare clean maps, we use both the MLR statistic and the 
NMSE (defined below). 

For an unpolarized and anisotropic gravitational wave 
background, the maximum-likelihood estimators of the 
signal-strength vector are given by 



k' 



(3.15) 



where S k are components of the dirty map (2.15) and 
V k are components of the deconvolved (clean) map. 
Note that the clean map, "P, are the values of V that 
maximize the statistic. We extend this single-baseline 
analysis to a multibaseline one by simply adding the 



dirty maps and beam matrices as = J2z=i $z an d 



®Nk> 



E N b r?k 
1=1 °1 



So, the maximum-likelihood estima- 



tors for a multibaseline network are given by 

V k = (.B*- 1 ) k k ,S%. 
We simulate the data with signal as [2] 



(3.16) 



xl(t,f)x 2 (t,f) 
(hl(t,f)h 2 (t,f)) 



(K(t,f)h2(t,f))+nUt,f)Mt,f), 
S fr HQf^V* 7(6*,*, l/l), (3.17) 



where V % is the injected source strength at the i th pixel, 
and 7 is the direction-dependent overlap reduction func- 
tion. We use a conjugate gradient (CG) method to solve 
the set of linear equations (3.15) and (3.16). 

The mismatch between two maps, injected and esti- 
mated, is measured using the normalized mean square 
error, 



NMSE := 



\v-v\ 



(3.18) 



12 




-0.2 -0.1 0.1 0.2 0.3 0.4 






-0.2 -0.1 0.1 0.2 



-0.2 -0.1 0.1 0.2 



(a) Injected map 



(b) H1L1 



(a) H1L1 



(b) L1V1 







-1 -0.5 0.5 1 1.5 2 



-0.2 -0.1 0.1 0.2 



-0.2 -0.1 0.1 0.2 



(c) L1V1 



(d) H1V1 



(c) H1V1 



(d) H1L1V1 



FIG. 10: Toy model of an extended source is shown in (a). 
Dirty maps made from simulated data containing signal from 
that source are shown in the last three panels for the three 
LIGO- Virgo baselines. 



FIG. 12: Difference between the clean maps of Fig. 11 and 
the injected map of Fig. 10(a). 
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(a) H1L1 (NMSE=0.4311) (b) L1V1 (NMSE=0.4171) 
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(c) H1V1 (NMSE=0.4735) (d) H1L1V1 (NMSE=0.2206) 

FIG. 11: Clean maps obtained by the deconvolution of the 
dirty maps of Fig. 10, using 20 CG iterations, are shown here. 



FIG. 13: The toy model of a localized source is shown in (a). 
Dirty maps made from simulated data from three LIGO- Virgo 
baselines are shown in the last three panels. 



1. Dirty maps and clean maps 

We performed numerical comparison of map-making 
performance for two types of toy sky patterns - (i) ex- 
tended, multideclination sky (Figs. 10 and 11), mimick- 
ing the (partially masked) image of the sky constructed 



by the WMAP satellite [ ], where essentially a modi- 
fied galactic structure stands out; and (ii) a relatively 
localized source peaked at the north pole (Figs. 13 and 
14). 

In both cases, the dirty maps from different base- 
lines are quite distorted compared to the injected maps. 
(Compare the last three plots with the first one in Fig. 
10 and in Fig. 13). However, the deconvolution proce- 
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FIG. 14: Clean maps obtained by the deconvolution of the 
dirty maps of Fig. 13, using 20 CG iterations, are shown here. 



individual baselines. 



2. Maximized-likelihood-ratio statistic 

We finally compute the MLR statistic introduced in 
Section II C as a figure of merit and, also, to demonstrate 
how this statistic can be powerful in identifying signal in 
noisy maps. The MLR statistic for both dirty and clean 
maps for the two types of sources considered here have 
been listed in Tables I and II. 

TABLE I: MLR statistic of dirty maps (A) versus clean maps 
(A c ) for the simulated maps in Figs. 10 and 11. 



Baseline 


A 


A c 


H1L1 


785.555 


783.271 


L1V1 


359.004 


358.940 


H1V1 


315.717 


315.662 


H1L1V1 


919.594 


917.600 
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FIG. 15: Difference between the clean maps of Fig. 14 and 
the injected map of Fig. 13(a). 



dure yields reasonably resolved maps for all the baselines 
(Figs. 11 and 14), signifying that none of the beam matri- 
ces are completely degenerate. The clean maps from the 
network are, however, of better quality, as can be seen 
from the corresponding low NMSE. To demonstrate this 
visually, we also show the difference between the clean 
and injected maps for the respective cases in Figs. 12 
and 15. As expected, the difference maps for the net- 
work look less noisy and uniform across the sky than the 



TABLE II: MLR statistic of dirty maps (A) versus clean maps 
(A c ) for the simulated maps in Figs. 13 and 14. 



Baseline 


A 


Ac 


H1L1 


284.652 


284.173 


L1V1 


39.308 


39.377 


H1V1 


64.129 


64.113 


H1L1V1 


294.419 


293.961 



It is intriguing to note that, for both dirty and clean 
maps, one obtains similar values. This suggests that a 
deconvolution effected with only a few tens of conjugate- 
gradient basis vectors does not cause a significant amount 
of information loss. 

To understand the significance of the MLR statistic 
in the present context, we perform two more exercises. 
First, we study the no-injection case; that is, we make 
dirty maps of simulated noise (Fig. 16), deconvolve it 
(Fig. 17), and obtain its MLR (Table III). The similarity 
of the values of this statistic for the dirty and clean maps 
proves the unitarity of our deconvolution method. One 

TABLE III: MLR statistic of dirty maps (A) versus clean maps 
(A c ) for simulated noise (which actually has an extremely 
weak signal added) in Figs. 16 and 17. 



Baseline 


A 


Ac 


H1L1 


0.512 


0.433 


L1V1 


-1.549 


-1.542 


H1V1 


1.105 


1.120 


H1L1V1 


0.208 


0.149 



can see that the MLR statistic (Table III) is small (~ 1) 
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(a) H1L1 (b) L1V1 




-0.1 -0.05 0.05 0.1 



(c) H1V1 

FIG. 16: Dirty maps made for simulated noise, without any 
injected signals. 




-0.06 -0.04 -0.02 0.02 0.04 0.06 -0.2 -0.1 0.1 0.2 

(a) H1L1 (b) L1V1 




-0.2 -0.1 0.1 0.2 -0.06 -0.04 -0.02 0.02 0.04 0.06 



(c) H1V1 (d) H1L1V1 

FIG. 17: Clean maps obtained by the deconvolution of the 
dirty maps of Fig. 16, using 20 CG iterations. 

in all these cases. We then introduce a small signal - 
the same as the extended (galaxylike) source considered 
before, but at a much reduced strength. 

Visually, the dirty maps are now weaker (Fig. 18), and 
clean maps almost do not show the obvious presence of 
any source (Fig. 19), but the MLR statistic (Table IV) 
provides a clear and reliable indication of the presence 
of a signal, thus proving its usefulness in the search for 
signal in a noisy map. 
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FIG. 18: The toy model of a very weak extended source is 
shown in (a). Dirty maps made from simulated data from 
three LIGO- Virgo baselines are shown in the last three panels. 
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FIG. 19: Clean maps obtained by the deconvolution of the 
dirty maps of Fig. 18, using 20 CG iterations. 

IV. DISCUSSION 

The search for an anisotropic stochastic gravitational 
wave background plays an important role in present GW 
research. In addition to setting interesting upper lim- 
its on astrophysical and cosmological backgrounds, the 
simplicity of the concomitant analysis reveals invaluable 
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TABLE IV: MLR statistic of dirty maps (A) versus clean maps 
(A c ) for the simulated maps in Figs. 18 and 19. 



Baseline 


A 


A c 


HlLl 


98.643 


98.670 


LlVl 


43.520 


43.536 


HlVl 


40.432 


40.475 


H1L1V1 


115.132 


115.176 



knowledge about the coherent performance of the GW 
detector network. 

So far, detailed analysis strategies have been developed 
to search for anisotropic background in pixel and spheri- 
cal harmonic spaces, and a general maximum-likelihood- 
based framework has been established to search in any 
convenient basis. The spherical harmonic search has been 
demonstrated using a network of detectors [ ]. In this 
paper, for the first time, we numerically implement the 
directed radiometer search, including deconvolution, for 
a network of detectors. These methods, in the past, were 
focused primarily on showing that the observed map is 
consistent with Gaussian noise or in estimating sky maps. 
The latter required the inversion of the convolution equa- 
tion, which itself assumed the network of detectors to be 
nondegenerate. Neither of these methods may work in 
the presence of excessive noise and weak signal. Most 
importantly, a statistically meaningful, all-sky combined 
statistic, in the form of an optimal "detection statistic," 
was needed in order to make precise statements about 
the presence or absence of a given background model in 
a map. Here, we proposed a MLR statistic, which yields 
a single number when computed on the dirty or the clean 
map and can be used as a detection statistic. By com- 
puting the MLR statistic for a couple of toy models of 
the background, we observe that the detection statistic 
is much larger than the noise-only case, even in the pres- 
ence of weak signals that are barely visible in dirty or 
clean maps. We corroborated these statements with re- 
sults obtained from extensive Monte Carlo simulations of 
a diffuse background of known shape in an ensemble of 
noise realizations. However, a more detailed study us- 
ing signals from a variety of background models is surely 
worth pursuing in order to determine how accurate the 
templates need to be in order to extract meaningful in- 
formation from weak backgrounds. 

We also compared the performance of individual base- 
lines and the whole network for the directed radiometer 
search using different figures of merit. Evaluating the 
performance of a network of GW detectors in SGWB 
searches is relatively straightforward compared to other 
GW signal searches [21-26]. This exercise was useful 
in drawing insights about the characteristics of a net- 
work that are particularly helpful in boosting its per- 
formance. Our overall observation, not surprisingly, is 
that the network improves performance in mainly three 
ways, namely, (1) by increasing the sensitivity by ob- 



serving each direction a greater number of times, (2) by 
observing the sky more uniformly, and (3) by probing 
each direction on the sky with additional detectors on 
the globe. The latter two enhancements lead to better 
localization of pointlike sources. This can be understood 
via the behavior of the Fisher information matrix: More 
detectors reduce its degeneracy and improve the well- 
posedness of the inverse problem. This, in turn, leads to 
a more accurate production of clean maps. 

Another question worth addressing in the future is 
about how closely spaced must the templates be on the 
parameter space to maximize the chances of detection 
with available computational resources. Indeed, the pro- 
posal for templated searches for SGWB signals is not new 
to this paper. For example, it has been addressed earlier 
in the context of isotropic searches (see Ref. [ ] and 
the references therein). Reference [ ] also introduced a 
metric on the parameter space of those signals so as to 
enable an experimenter to infer what the principle axes 
are on that space and how fine a template bank one can 
afford based on the computational resources available. A 
similar study can be carried out for finding a more opti- 
mal spacing of templates for directed searches than the 
one used here. 

Whereas results presented here were derived for Gaus- 
sian noise, the codes used can be applied to real data as 
well. Indeed, the performance of the proposed statistic in 
real data sets from the LIGO and Virgo detectors can be 
determined through hardware injections that were done 
in the recent science runs, such as the ones described in 
Ref. [ ], and supplementing them with multiple soft- 
ware injections to improve the statistics. The expected 
improvement of network sensitivity over individual base- 
lines, as demonstrated here, merits the investment re- 
quired for extending the current single-baseline analysis 
efforts [ , ] to a multibaseline one. This conclusion is 
strengthened by the fact that adding a detector to a base- 
line can potentially mitigate the contribution of cross cor- 
related environmental noise that affects only one of the 
three resulting baselines. Including VI, which is on a dif- 
ferent continental plate than the HlLl baseline, can serve 
this purpose. Employing a null-stream statistic [26, 27] 
to complement the detection statistic might also help in 
discriminating against such noise. 



Acknowledgments 

We thank Joe Romano, Stefan Ballmer, and War- 
ren Anderson for discussions, careful reading of the 
manuscript, and helpful comments. We would also like 
to thank Bruce Allen, Sanjeev Dhurandhar, Albert Laz- 
zarini, Vuk Mandic, Himan Mukhopadhyay, Alan We- 
instein, and Holger Pletsch for helpful discussions. SM 
would like to acknowledge the Centre National d'Etudes 
Spatiales (France) for supporting part of the research. 
This work was supported in part by NSF Grant No. 
PHY-0855679. Part of the research described in this 



16 



paper was carried out at the Jet Propulsion Labora- 
tory, California Institute of Technology, under a contract 
with the National Aeronautics and Space Administra- 
tion. LIGO was constructed by the California Institute 
of Technology and Massachusetts Institute of Technology 
with funding from the National Science Foundation and 
operates under cooperative agreement PHY-0757058. 



2 ( B (fc)(fc)) [ 
d 



Bk'k 



k'=k 



d 2 



d\± k >d(\) k 
d 



k'k 



k'=k 



k'=k 



(A5) 



Appendix A: Parameter accuracy 

The match can be rewritten as 
M = 1 
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The components of the above g a $ matrix are obtained 
from the derivatives of the beam matrix: 
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The estimation error (3.11) is obtained from Eqs. (A2) 
and (3.10) by utilizing the fact that the inverse of the 
determinant of a matrix is the same as the determinant 
of the inverse of that matrix. 
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